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First we briefly review our covariant Hamiltonian approach to quasi-local energy, noting 
that the Hamiltonian-boundary-term quasi-local energy expressions depend on the cho- 
sen boundary conditions and reference configuration. Then we present the quasi-local 
energy values resulting from the formalism applied to homogeneous Bianchi cosmologies. 
Finally we consider the quasi-local energies of the FRW cosmologies. Our results do not 
agree with certain widely accepted quasi-local criteria. 
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1. Introduction 



The localization of energy-momentum for any gravitating system (and thus for all 
physical systems) is still an outstanding fundamental problem. In view of conserva- 
tion, and the fact that sources exchange energy-momentum locally with the grav- 
itational field, some kind of local description for gravitational energy-momentum 
was expected. However all attempts at constructing such an expression led only 
to reference frame dependent quantities, generally referred to as pseudotensors.'^'^ 
It became apparent that the gravitational field itself, unlike all matter and other 
interaction fields, has no proper energy-momentum density. This fact can be under- 
stood as a consequence of Einstein's equivalence principle. 51 The energy-momentum 
of gravity — and thus the energy-momentum for all physical systems — is inherently 
non-local. The modern idea is quasi-local: energy-momentum is associated with a 
closed surface bounding a region. 
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Many quasi-local expressions have been proposed, but presently there is no con- 
sensus as to which is the most suitable, or even as to which properties a good 
expression should have.Sl Various lists of desiderata for a ''physical" quasi-local en- 
ergy have been presented; according to a well-known one—' the quasi-local energy 
should be 

• zero for flat space, 

• for spherical symmetric ~ standard value, 

• the ADM mass at spatial infinity, 

• the Bondi mass at null infinity, 

• for the apparent horizon ~ standard value, 

• positive. 

Our Hamiltonian based quasi-local results do not satisfy the first and last of these 
criteria. There is a stronger forni'^ of the first requirement, namely that the energy 
vanish iff the quasilocal region is flat space. Our analysis of the quasi-local energy 
of cosmological regions leads us to propose a certain modification of this stronger 
form. 



2. The covariant Hamiltonian approach 

Here we briefly sumn iarize the re levant parts of our covariant Hamiltonian approach 
to quasi-local energy.l2 E l 7 | 8 | 9 | 10 | 



2.1. First order Lagrangian 

The first order Lagrangian for an f-form field and its conjugate momentum p has 
the form 

C = d(p hp - K{if,p). (1) 

The variation of this 4-form, 

6C = d{5ip Ap)+Siph — + — A6p, (2) 
dip op 

leads to the first order equations of motion 
X r r\ r 

— := - 9pA = 0, —~ -^dp - d^h = 0, (3) 
op dip 

where ? := (—1)-'^. 

As a simple example of this formalism consider electromagnetism. The first order 
Lagrangian 4-form for the (source free) U(I) gauge field one-form A is 

Cem = dAAH -^*H ah. (4) 
Variation leads to the pair of first order equations 

dH ^0, dA-*H = 0. (5) 
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From F :— dA — *H one finds H — — *F, hence the first equation becomes d*F = 0. 
Thus we have the vacuum Maxwell equations. 



2.2. Translation invariance and Noether current 

Infinitesimal diffeomorphism invariance (in terms of the Lie derivative) requires that 
(|2l) becomes an identity under the replacement S —^ £n' 

S jO S 

diisiC = £ njC = d(£ Nip a p) + £nV A — h — A £nP- (6) 

oiy9 op 

This simply means that £ is a 4-form which depends on position only through the 
fields ip, p. For this to be the case the set of fields in C necessarily includes dynamic 
geometric variables, which means gravity. 

From this identity it follows that the "translational current" density (3-form) 

n{N) ^ £NipAp-iNC (7) 

satisfies the "conservation law" 

^dH{N) = £n(pA^ + ^ /\£nP- (8) 
0(p op 

Consequently, "on shell" (i.e. when the field equations are satisfied), the integral of 
the current over a spatial region will give a conserved quantity for each vector field 
N. Note that, just like other Noether conserved currents, H{N) is not unique: it 
can be modified by adding the differential of any 2-form. 

With geometric gravity included, we have also local diffeomorphism invariance, 
which gives rise (in accordance with Noether's second theorem) to a differential 
identity. Explicit calculation shows that 7i{N) ~ £n^ A p — always has the 
form 

n{N) ^ Nt'n^ + dBiN). (9) 

Thus we find that diNf^H^, + dB{N)) = dN^' A + TV^dH^ is proportional to the 
field equations, therefore Hfj. vanishes "on shell" . Hence for gravitating systems the 
Noether translational "charge" — energy- momentum — is quasi-local: it is given by 
the integral of the boundary term, B{N). But this boundary term as noted can be 
completely modified to any value. The Hamiltonian approach includes an additional 
principle which naturally tames this ambiguity.121 



2.3. Hamiltonian approach 

Energy can be identified as the value of the Hamiltonian associated with a time- 
like displacement vector field TV. The Hamiltonian H{N) is given by an integral 
of a suitable Hamiltonian 3-form (density) 7i{N) over a 3-dimensional (space-like) 
region E. Generalizing L = qp — H , from the first order Lagrangian one constructs 
the Hamiltonian 3-form by projecting along a "time-like" displacement vector field: 

InC^ £n^ ^p-niN). (10) 
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The Hamiltonian density thus turns out to be just the Noether translational cur- 
rent (O identified above. As aheady noted it satisfies the relation ([9]) with Hfj. 
vanishing on shell. Consequently the quasi-local energy — regarded as the value of 
the Hamiltonian — is then determined only by the boundary integral: 



The two parts of the Hamiltonian have distinct roles. The 3-form, 7i^, although it 
has vanishing numerical value, generates the equations of motion. For our concerns 
here the Hamiltonian boundary term B{N) is the key quantity. It plays a dual role: 
determining both the the quasi- local values and the boundary conditions. 

2.4. Quasi-local quantities 

The Hamiltonian boundary term B{N) determines the various quasi-local values 
corresponding to the Poincare transformation of space-time: 

• Energy < — > N a time- like displacement, 

• Linear momentum < — > a spatial translation, 

• Angular momentum < — > a rotation, 

• center-of-mass moment •* — > a boost. 

However we noted that B{N) can be adjusted; then it would give different conserved 
values. What do they all these different values mean physically? 

2.5. Boundary Conditions 

The variational principle contains an additional (largely overlooked) feature which 
distinguishes all of these choices: the boundary variation principle, i.e. the boundary 
term in the variation tells us what to hold fixed on the boundary — it determines the 
boundary conditions. Different Hamiltonian boundary term choices are each associ- 
ated with distinct boundary conditions. (In this way this formalism gives a s pecifi c 
physical significance to each of the traditional energy-momentum complexes .^^^^^ 




This feature is similar to that of some familiar physical systems. For example 
in thermodynamics the suitable measure of energy: the internal energy, enthalpy, 
Helmholtz, or Gibbs free energy depends on the system's boundary conditions. An- 
other good example concerns moving a dielectric within a parallel plate capacitor. 
The work needed, and thus the appropriate energy density expression (the sym- 
metric or the canonical tensor) depends on the boundary condition: fixed charge 
or fixed potential. Thus one can see that there always are various distinct physical 
"energies" which correspond to how a system interacts with the outside through its 
boundary. 




(11) 




I 



I 



1, 2008 6:3 WSPC/INSTRUCTION FILE CosPA06 



Quasi-local energy for cosmological models 5 

2.6. Reference Configuration 

In general (in particular for gravity) it is necessary (technically, in order to guarantee 
functional differentiability of the Hamiltonian on the phase space with the desired 
boundary conditions) to adjust the boundary term, B{N) — i^ip A p, which is 
naturally inherited from the Lagrangian. Moreover a reference configuration, (p and 
(which determines the ground state) is essential (especially for gravity where the 
ground state is not vanishing field but rather the Minkowski metric) in particular 
to allow for the desired phase space asyniptotics. 

2.7. Quasi-local Expressions 

With Aip :~ ip — (fi, Ap :— p — p, where the bar indicates the reference value, we 
found two boundary choices (essentially Dirichlet and Neumann) which have the 
indicated covariant boundary terms in STi.. 

By, ^ iNf A Ap - <;Alp AiNp, ~^ iN{S(pAAp), (12) 
Bp — iNif A Ap — i^Aip A inP, —iN{AipA6p). (13) 

3. Application to GR 

For Einstein's (vacuum) gravity theory, General Relativity (GR) a first order La- 
grangian is 

^GR = T^i?"/3 A ?7„^ (14) 

where F"^ is the connection one-form, R"^ := dT^/j + T"^ A is the curvature 
2-form and 77"^ := *(i9" A t?^). 

Our general formalism with — > F"^ and p ^ rj^^ gives the quasi-local expres- 
sions for GR. We have two e xpre ssions for different types of boundary conditions. 
One of the choices stands out.l^ 

B{N) ^ [AF"^ A z^t?/ + {DN)p^Ar^J] . (15) 

This is a Dirichlet type condition for a covariant object, the orthonormal frame 
field. Asymptotically this expression gives not only the ADM (spatial infinity) and 
Bondi energy (null infinty) but also the Bondi energy flux. Moreover this expression 
is distinguished by satisfying a positive energy property. 

In the cases considered here, the contribution of the second term in (jlSp vanishes. 

4. Homogenous Cosmologies 

Homog ene ous cosmologies (non-isotropic in general) are described by the Bianchi 
models!^ the orthonormal coframe has the form 

d°^dt, 19'^ = /^%(^)^7^ (16) 
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where the spatially homogeneous frames satisfy 

da*^ = ic^jcr* A(T^ (17) 

where the C'^ij are certain constants. The associated space-time metric is thus 

ds2 ^ _^f2 ^ g^.(t)a'(^x)a^{x), (18) 

where gij := 5abh°'ih^j (which need not be diagonal). 

There are 9 Bianchi types distinguished by the particular form of the structure 
constants C'^y , especially by the value oi Ak := C^ki- They fall into two special 
classes: 

• class A {Ak = 0): Types I, II, VIo, VIIo, VIII, IX; 

• class B {Ak ^ 0): Types III, IV, V, VI/,., VII,,. 

The respective scalar curvatures are: vanishing for Type I, positive for Type IX, 
negative for all the other types. It should be mentioned that certain special cases 
can be isotropic, specifically isotropic Bianchi I, V, IX are equivalent to the usual 
FRW k = 0,-1, +1. 

For the natural choice of = 9*, the Dirichlet type boundary condition, the 
Bianchi homogenous frame as boundary value, and with the reference being the 
static homogenous cartesian frame, the energy within a spatial volume V according 
to our favored quasi-local expression igl^ 

E{V)^j^A,Akg^\t)V{t)>Q. (19) 

The result is true for all regions and for all types of sources including dark matter, 
dark energy a/o a cosmological constant. More specifically it vanishes for all class 
A models and is positive for all class B models. Note: this is entirely consistent with 
the important requirement that E = Q for closed universes, since all homogeneous 
class A models can be compactified and class B models cannot .^I^ 

5. FRW cosmology 

The Friedman- Robertson- Walker (FRW) (homogeneous and isotropic) metrics have 
the form 

ds^ = -dt^ + a^{t)dl'^. (20) 

The spatial metric dP has constant curvature. The FRW spatial metric has several 
equivalent manifestly isotropic-about-a-chosen-point forms: 

di^ = dp^ + Y.'^dn^ = + r^dn"^ = ^ (dR^ + R^dn^) , (21) 

1 - kr^ [1 + (A:/4)i?2]2 ^ / ' ^ ; 

where S = (sinhp, p, sinp) for k = (— 1,0, +1), respectively. 

A natural choice in this case isN = dt, Dirichlet type boundary conditions, FRW 
frame boundary values, with the reference being the static flat cartesian frame. The 
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energies within a fixed radius for the three FRW cases can be represented in several 
equivalent forms (their identity follows from S = r = -R/(1 + kFt? /A)): 

akR^ 



E = aS(l - E') = or 1 - ^1 - kr'^ 
More specifically, 
£'fe=_i = asinhp(l — coshp) = ar 
Ek=o = 0, 



Ek=+i = asinp(l — cosp) = ar 



1 - 



[l + (A;/4)i?2]2- 

~ 2(l-i?2/4)2 
2(l + i?2/4)2 



(22) 



<0, 



> 0. (23) 



6. Discussion 

According to our favored quasi-local energy expression, homogeneous choices give 

vanishing energy for all regions of Bianchi class A models and positive energy for 
class B. Isotropic choices give energies proportional to the spatial curvature param- 
eter k: vanishing for flat, negative for the open model, and positive for the closed 
model (but nevertheless vanishing, as required, when the considered volume is ex- 
tended to include the whole universe). 

Some of the Bianchi models can be isotropic, specifically 

• isotropic Bianchi I (class A) = FRW/c=o; 

• isotropic Bianchi IX (class A) = FRWfe=+i; 

• isotropic Bianchi V & Yllh (class B) = FRWfe=_i. 

Note that our quasi-local expression thus can give different energy values to exactly 
the same geometry. This is not at all mysterious; it is clearly a consequence of dif- 
ferent reference and boundary value choices. Homogenous boundary values are not 
the same as isotropic boundary values. To understand the physical and geometric 
meaning of the differences between the homogeneous and isotropic choices in detail, 
we need to do more calculations using the rather complicated relations between the 
FRW and Bianchi coordinates. Meanwhile from our analysis it seems that the homo- 
geneous choice is more suitable physically than the isotropic-about-a-chosen-point 
choice in general, since it gives a non-negative energy. 

It is also noteworthy that for the case of the open FRW (fc = —1) with vanishing 
matter, the solution to the Einstein equation gives a{t) = t. It can be directly verified 
that the geometry is then really just Minkowski space in non-standard coordinates, 
yet our expression gives a non-vanishing energy, which, moreover is negative with 
the FRW choices. 

Concerning two of the quasi-local desiderata, for the expression and bound- 
ary/reference choices considered we found that 

• positivity need not hold; 

• "zero energy iff flat Minkowski space" need not hold in either direction. 
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It seems that these quasi-local criteria should be reconsidered. 

As we prefer positivity, we are inclined to see our negative result as disfavoring 
the FRW isotropic-about-a-chosen-point boundary/reference choice. 

Regarding E = 0, clearly any homogeneous measure of quasi-local energy in 
Bianchi I models must necessarily vanish for all regions — since these models can be 
compactified (with 3 torus topology identifications) on any scale, and the energy of 
a closed universe must vanish. In light of this, we propose that the "unique quasi- 
local E = Minkowski ground state" requirement be replaced by something like 
''E{V) = iff a neighborhood of V can be compactified." 
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